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The successful production of Bose-Einstein conden- 
sates (BEC) has led to measurements indicating that 
BECs have coherence properties similar to lasers. This 
suggests that the next stage in obtaining precision mea- 
surements in atom optics is the production of atom beams 
with sub-Poissonian atom statistics - as has been widely 
demonstrated in photonics applications. Indeed, there is 
already some indirect evidence of atom-number squeez- 
ing, but experiments and theoretical proposals to date 
have focussed on trapped condensates and phase sensitive 
measurements [|l|, ^, ^. A possible route towards more 
robust phase- insensitive applications of atomic squeezing 
is to produce quantum correlated atom laser beams in 
which the correlations are directly built in the statistics of 
the particle numbers, as in parametric down-conversion 
in quantum optics Q . Quantum correlated or entangled 
photon pairs from parametric down-conversion are one of 
the most powerful resources of quantum optics, and our 
motivation here originates from the intriguing prospect of 
possessing a matter-wave analog of this resource. As well 
as the immediate possibility of improved atomic interfer- 
ometry, an exciting prospect would be the development 
of new tests of quantum measurement theory for massive 
particles with space-like separations, since all previous 
tests using down-conversion methods were restricted to 
massless photons. 

In this paper we propose a robust scheme for achiev- 
ing strong quantum correlations between two counter- 
propagating atomic beams, relying on the process of 
photo-dissociation of a molecular Bose-Einstein conden- 
sate. Experiments towards production of molecular con- 
densates H, ^], together with a number of theoretical 
studies of coupled atomic-molecular BEC systems are 
the subject of much intense activity at present 0, I). 
We anticipate that the formation of a molecular BEC is 
a matter of time, and consider it as the starting point 
for a twin atomic beam experiment with relative atom 
number fluctuations reduced below the level predicted by 
either thermal or coherent (Poissonian) statistics. The 
method automatically yields two counter-propagating 
beams through momentum conservation, and is robust 
against changes in mode structure, coupling constants, 
or even absorptive losses, provided they are small. 

An important feature of our proposal is that it does 
not rely on atomic intcrfcromctry or local oscillators to 



generate the resulting correlations, which makes it more 
practical than recent related proposals [§, ^ . In addition, 
the present scheme is not susceptible to phase noise from 
self-phase modulation. In the following analysis, we first 
analyze a simplified theory of a uniform, non-depleted 
molecular condensate, then include the effects of molecu- 
lar trapping, depletion, phase diffusion and atomic losses. 

The quantum field theory effective Hamiltonian for the 
atomic {^i) and molecular {'^2) fields, taken for simplic- 
ity to be confined to one space dimension, is: 

with the commutatotion relation [^'^(a;, i), i)] = 

5ij6(x~x'). Here stands for the usual kinetic energy 
term, Vi{x) is the trap potential (including internal en- 
ergies), Uii ~ Anhai / (Ami) is the atom-atom coupling 
constant in one dimension, where mi is the mass, a is the 
three-dimensional 5'-wave scattering length, and A is the 
confinement area in the transverse direction, with similar 
results for the molecule-molecule and molecule-atom in- 
teractions. The term proportional to x{t) describes a co- 
herent process of molecule- atom conversion via either one 
or two-photon (Raman) photo-dissociation, where x{t) is 
the coupling constant related to the transition matrix el- 
ement(s) and the amplitude(s) of the photo-dissociation 
laser (s) which have an overall frequency (difference) uj. 

The time dependence of x{t) controls the duration of 
dissociation process and will be set to x{t) = Xod{t)8{ti — 
t), where is a step function. This means that the dis- 
sociation can be stopped after a short time interval ti 
followed by free evolution of the atomic field. Starting 
from a pure molecular condensate, the molecular field 
can be initially described semiclassically, via its initial 
coherent amplitude '^2{x,0) = {'^2{x,0)). We assume 
'n2{x) — \'i>2{x,0)\'^ is the initial molecular BEC density 
in a harmonic trap in the Thomas- Fermi limit, with the 
axial half-length xq- 

The atoms are assumed to be untrapped longitudinally 
(they may be in an m = magnetic sublevel) yet confined 
transversely (they may be in a transverse optical trap). 
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so that the atomic field can effectively be treated as a free 
one-dimensional field, initially in a vacuum state. One- 
dimensional trapping of condensates has been achieved 
experimentally |lO[ , so this is not unrealistic. In what fol- 
lows, we will neglect the atom-atom collisions for simplic- 
ity, owing to the fact that we restrict ourselves to short in- 
teraction times during which the atomic density does not 
grow to high values, so that the self-interaction term pro- 
portional to Uii can be neglected. In addition, we choose 
the absolute value of the detuning | A| to be always non- 
zero and much larger than the magnitude of C/i2(^J^2) 
so that the atom-molecule S -wave scattering is negligi- 
ble too. The detuning A = Vi{0) - [V2{xo) + uj]/2, where 
14(a^o) = ^2(0) -I- 1/22'^2(0), is proportional to the energy 
mismatch between the atomic and molecular fields. 

The system has a direct analogy with travelling-wave 
parametric down-conversion in non- linear optics 
Here the role of xo is played by a non-linear crystal with 
second-order susceptibility x^'^\ ^ind a finite interaction 
time is analogous to the crystal length L = vtiin the di- 
rection of propagation, where v is the group velocity of 
the fundamental (higher frequency) beam. The detuning 
A is analogous to the optical phase mismatch, while the 
atomic kinetic energy is analogous to dispersion. 

To proceed with the analysis we introduce a character- 
istic length scale do and time scale to = 2mid^/h. Next, 
transform to dimensionless fields, in rotating frames: 

MLr) = 7^:^*2(2;, i)e*^^(^«)S (2) 

where ^ = x/cIq and r — t/to are the dimensionless co- 
ordinate and time. We also introduce dimensionless de- 
tuning 6 — Ato and coupling K{t) = x{t)to/Vdo- 

The Heisenberg equations of motion for the field oper- 
ators, in dimensionless form are: 



dr 
Ot 



' de 



(3) 



We have introduced an effective potential ^2(0 — 
[V2{£.do) - V2(Corfo)]^o + w'/'2V'2, where u = U22to/dQ, for 
notational simplicity. 

To gain some insight into the underlying physics of 
correlated atomic beams, we first consider an idealized 
and analytically solvable model corresponding to an un- 
depleted and uniform molecular condensate at density 
712(0), that fills the entire space from —1/2 to 1/2, with 
periodic boundary conditions at —1/2 and 1/2. The atom- 
molecule coupling X = Xo is assumed to be constant dur- 
ing the whole evolution time from to r. In this case 
the coherent amplitude of the molecular field can be ab- 
sorbed into an effective gain constant g = Ko-y/n2(0)(io 



(where kq = Xoto/Vdo) which we assume without loss of 
generality to be real and positive. 

Solutions to the resulting linear set equations of motion 
for the atomic field are easily found in momentum space, 
where we expand ipi{^,T) in terms of single- mode anni- 
hilation operators: ?/'i(^,r) = J^q^qi''')^^'^^ where 
q — do k is a dimensionless momentum. The single- 
mode operators aq satisfy the usual commutation rela- 
tions [aq{T),al^,{T)] = ^qq'- The corresponding Heisen- 
berg equations of motion have the following solution 

dqir) = Aq{T)aq{0)+Bq{T)alqiO), 

alq{T) = Bq{r)aq{0) + A;{T)alq{0). (4) 

where 



Aq{T) ^ cosh (gqT) - iXqSinh (ggT) /gq , 
Bqir) = g smh {gqT) / g^ , 



(5) 



with \q = q 



(5, and gq = (g^ — A^)^/^. Solutions 
of this type to the classical counterpart of the operator 
equations are well known in optics |pl| , while in quantum 
optics the operator equations in the context of squeezing 
of nonlinear propagating fields were studied in . 

Knowledge of the initial state of the atomic field, which 
is the vacuum state with 6,(0) |0) = 0, allows us to cal- 
culate any operator moments at time t. The parameter 
gq is the gain coefficient; if real, it causes a growing cor- 
related output for the momentum component q, while if 
imaginary it leads to oscillations. 

For example, the particle number distribution in mo- 
mentum space is given by {a)q{T)dq(T)) = Bq{T). For 
S < 0, the function has two distinct global max- 
ima located at g-values where Xq — 0. This gives the 
two most probable momentum values 50 = ± 'v/RI? corre- 
sponding to a zero effective phase mismatch term (q'^+S), 
provided 6 < 0. The total average number of atoms 
{N{t)) = J2q {dl{T)aq{T)) is glvcu by 



{N{r))=J2 {9/gqfsmh\gqT), 



(6) 



which grows exponentially with r. 

To analyze correlations and relative number squeezing, 
we define particle number operators iV_(r) and N^{t) 
containing only negative or positive momentum compo- 
nents, respectively: iV_(+)(r) = J2q<o{q>o) 
We next consider the normalized variance V{t) of the 
particle number difference [N^{t) — iV_|_(T)], which - in 
normally ordered form - is given by: 

Vir) = 1 + (: [A(7V_ - N^)]^ :>/((^-) + (^+» (7) 

where AX = X — {X), and V{t) < 1 implies squeezing 
of fiuctuations below the coherent level - which is due to 
quantum correlations between 7V_ and N^. 
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Calculating (: (iV±)^ :) and {N^N+), and using the 



fact that B^-\Aq\ 



-1, gives 



V{t) = 1 + 



g>0 



q>0 



Bl^Q, (8) 



implying perfect (100%) squeezing in the particle number 
difference, at least in this idealized calculation. 

We now turn to the analysis of more realistic non- 
uniform case and including the molecular field depletion, 
as described by the original coupled operator equations 
for the fields, Eqs. (|). These are solved numerically 
using equivalent stochastic (c-number) differential equa- 
tions in the positive-P representation , where we ad- 
ditionally include a coupling to an atomic loss reservoir 
to describe possible linear losses at a rate 7: 



dip! 
dr 

dr 

dr 

dr 



i d^ip2 



l2„/,+ 



Here V'i and are complex stochastic fields correspond- 
ing respectively to the operators ipi and ipj, v{^,t) — 
[V2{^do) — V2{£,odo)]to + u'4>2'4^2 represents the effective 
potential, and rji , rj'^ are four real independent delta- 
correlated Gaussian noise terms: ''')'7i '''')) ^ 

We consider molecules as an initial coherent field cor- 
responding to Thomas-Fermi inverted parabola for the 
molecular density. In this case, assuming that "02 (■^,0) 
is real, we have V'2(C,0) = a/ rfo'T-2(C)' where n2(x) = 
"2(0) [1 - (^/?o)^] 6'(Co - The molecular conden- 

sate axial half-length is denoted via ^0 which is deter- 
mined by the trap geometry, and we assume repulsive 
molecule-molecule interactions. The time duration for 
the molecule-atom conversion is controlled via k(t) — 
ko9{ti — t), so that k(t) = for t > ti. Once the 
dissociation is stopped, we continue the evolution of the 
resulting atomic field in free space to allow spatial sepa- 
ration of the modes with positive and negative q values. 

For spatially separated components, we can introduce 
a pair of particle number operators: 

N-i+){T) = r^^^ vii(c,r)^i(e,T)de (10) 

J-i/2(0) 

Next, we define the normalized variance V{t) of the par- 
ticle number difference [N-{t) — N+{t)] as before [see 
Eq. (^] , and evaluate the average values numerically us- 
ing the positive-P Eqs. (||) and the standard correspon- 
dence between the normally ordered operator averages 
and the c-number stochastic averages p3|. 
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FIG. 1: Atomic density ni(^,r), total number of atoms 
Ni{t), and the variance V{t) , for 40000 trajectory averages 
and for: S = -5.017 x lO" , u = 1.8, kq = 84, and 7 = 5 
. The dimensionless dissociation time in this simulation is 
n = 8 X 10"''. With the choice of the length scale do ~ 30 
/im and the mass of a *^Rb atom, we obtain the time scale 
to ~ 2.446 s, so that r = 0.011 corresponds to the total time 
window t ~ 27 ms, while the dissociation time ri scales to 
ti ~ 2 ms. 

The results for V{t), the total number of atoms in the 
two beams Ni{t) = {N^{t)) + (iV+(r)), and the density 
distribution ni(^,r) — (V'j (^, t)V'i (^, r)) are represented 
in Fig. ^ for parameters that are reasonable in current 
experiments and including a linear atomic loss term as 
discussed below. 

The initial growth of the fluctuations in the atom num- 
ber difference, during the time interval when the atom- 
molecule coupling is switched on, is due to the fact that 
the quantities {N-{t)) and (iV_|_(r)) each include atoms 
traveling in opposite directions inside the molecular BEC 
region. The density distribution in coordinate space 
at this stage is single peaked, due to the fact that it 
contains amplified contributions from both momentum 
components. The fundamental correlation of opposite 
momentum components is therefore not visible in the 
atomic density initially. Once, however, the interaction is 
switched off and the correlated atom pairs fly apart with- 
out further parametric amplification, there is a double- 
peaked distribution, and we see a rapid reduction of the 
variance below the coherent level, V{t) < 1. 

The physical reason for the correlation is momentum 
conservation, which requires that each emitted atom with 
q > be accompanied by a partner atom having q < 0. 
In order to conserve energy, this pairing only occurs for 
(5 < 0, which allows potential energy in the molecule to 
be converted to atomic kinetic energy for selected modes 
with q values around qo = ±^/\S\. 

As the scheme relies on conservation laws for its opera- 
tion, it should be insensitive to the exact mode-structure. 
Clearly, quasi-one dimensional traps are preferable for 
reasons of directionality, but we expect similar results 
even if there is no transverse trap, provided that the 
molecules are confined in a high aspect-ratio (cigar- 
shaped) trap to allow gain-guiding of the atoms. 
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Investigations of losses show that these have a mini- 
mal effect on squeezing, provided the atoms that are lost 
are only a fraction of the atoms produced. For exam- 
ple, including an atomic loss term at a rate 7/2 — 10, 
corresponding to the loss of ~ 10% of atoms during the 
free evolution time interval, gives V ~ 0.07 at r = 0.011 
as shown in Fig. 1. This corresponds to a rather high 
93%) degree of squeezing below the coherent level. 
Additional effects may occur when there is an atomic 
potential, or when there is a strong atom-molecule scat- 
tering, but this has very little effect on squeezing when 
\S\ 3> \ua\, if Wq is the scaled effective atomic potential. 

The parameter values of Fig. 1 are derived using 
a length scale of do = 30 /Ltm taken to be equal to 
the molecular condensate axial half-length R^, so that 
Rx = xq — do and ^0 = 1- We assume that the molecular 
BEC is formed in a highly elongated trap and at den- 
sities that satisfy the conditions for the crossover from 
3D to ID (see, e.g., @). Assuming that the molecule- 
molecule scattering length 02 is of the same order of mag- 
nitude as the scattering length of ®^Rb atoms, we take 
02 = oi = 5.4 nm. We also take the molecular BEC lin- 
ear density at the trap center 712(0) = 3.7 x 10^ m~^ and 
an aspect-ratio of 100. This implies that the condensate 
transverse radius is R± = 0.3 /xm, so that the 3D peak 
density is about 3.27 x 10^^ m~^. The initial total num- 
ber of molecules is N2 — 1.48 x 10"^. Using the value of 
R± to scale out the transverse confinement {A Tri?^), 
we can next estimate the one dimensional values of Uij 
and xo- The ID value of xo is obtained according to 
Xo = x'^^^yV'^Rl, where we take x^^^) = 2 x 10^"^ 
m^/^/ s [0. The molecular condensate half-length of 30 
fi m corresponds to the trap axial oscillation frequency 
ujx/2Tr = 4.1 Hz, where lOx = \/'2hU22n2{0)/{m2xl). We 
choose A — —2.051 x 10"* s~^, thus assuring that the 
the atom-molecule 5'- wave scattering is negligible in this 
calculation since |A| ^ |f^i2|(^2^2)i where U12 is esti- 
mated using ai2 ~ —9.25 nm |5|. The final values of the 
parameters 6, u, and kq that we arrive at are specified in 
the figure caption, together with the relevant time scales. 

In summary, we have shown that photo-dissociation of 
a molecular BEC can provide a simple yet robust scheme 
for quantum squeezing of relative number fluctuations in 
two counter- propagating atomic beams. The effects of 
molecular condensate trapping and depletion, molecular 
self-phase modulation, and atomic absorption have all 
been included in our calculations. Our method does re- 
quire high efficiency atom counting techniques Jl^ , which 
are currently the subject of intensive activity. 

Applications may emerge from the use of these 
quantum-entangled twin beams to produce a single beam 
with a well defined particle number, which can be 
achieved by a, destructive mea,siirement on the partner 



subtle experiments may be feasible in future, including 
possible demonstrations of Einstein-Podolsky- Rosen cor- 



relations or Bell inequalities |15 in matter-wave quadra- 
tures. Such experiments would open the way to novel 
tests of quantum mechanics for macroscopic numbers of 
massive particles. 
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useful discussions. 
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